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We theoretically study the competition between two possible exotic superconducting orders that 
may occur in graphene-like systems, assuming dominant nearest-neighbor attraction: the gapless 
hidden superconducting order, which renormalizes the Fermi velocity, and the Kekule order, which 
opens a superconducting gap. We perform an analysis within the mean-field theory for Dirac 
electrons, at finite-temperature and finite chemical potential, as well as at half filling and zero- 
temperature, first excluding the possibility of the coexistence of the two orders. In that case, we 
find the dependence of the critical (more precisely, crossover) temperature and the critical interaction 
on the chemical potential. As a result of this analysis, we find that the Kekule order is preferred 
over the hidden order at both finite temperature and finite chemical potential. However, when 
the coexistence of the two superconducting orders is allowed, using the coupled mean-field gap 
equations, we find that above a critical value of the attractive interaction a mixed phase sets in, in 
which these orders coexist. We show that the critical value of the interaction for this transition is 
greater than the critical coupling for the hidden superconducting state in the absence of the Kekule 
order, implying that there is a region in the phase diagram where the Kekule order is favored as 
a result of the competition with the hidden superconducting order. The hidden superconducting 
order, however, eventually sets in and coexists with the Kekule state. According to our mean-field 
calculations, the transition from the Kekule to the mixed phase is of the second order, but it may 
become first order when fluctuations are considered. Finally, we investigate whether these phases 
could be possible in honeycomb superlattices of self-assembled semiconducting nanocrystals, which 
have been recently experimentally realized with CdSe and PbSe. 

PACS numbers: 71.10.Li, 74.25.Dw 


I. INTRODUCTION 


Ever since its isolation, graphene and graphene-related 
topics have attracted much attention from both theoret¬ 
ical and experimental condensed-matter communities 
The monolayer of carbon atoms is light, transparent, 
flexible, strong and conductive, which makes it a per¬ 
fect candidate for industrial applications. The appear¬ 
ance of these properties in the material are for a large 
part due to the assembly of the atoms in a honeycomb 
geometry, which consists of two interpenetrating trian¬ 
gular Bravais lattices, referred to as sublattices A and 
B. The lattice symmetry together with time-reversal 
gives rise to the hallmark feature of the graphene sys¬ 
tem - pseudo-relativistic massless Dirac fermions, which 
are low-energy quasiparticles close to the Dirac points lo¬ 
cated at the two inequivalent, time-reversal conjugate K 
and K' = — K momenta at the corners of the Brillouin 
zone (BZ). As a consequence, the density of states lin¬ 
early vanishes close to the Dirac points. The semimetal- 
lic ground state is therefore protected against the effects 
of weak interactions, and an intrinsic superconducting 
state in half-filed graphene could thus only be possible for 
sufficiently strong attractive interactions.!^ Onsite attrac¬ 
tion supports an s—wave spin-singlet superconducting 
stat e,!^ whereas nearest-neighbor attraction may lead to 


the formation of the Kekule state, which breaks transla¬ 
tional lattice symmetry and opens up a gap at the Dirac 
points,!^ while, at the same time, a gapless superconduc¬ 
tor may also set inW Therefore, it is of fundamental im¬ 
portance to address the competition of these two super¬ 
conducting orders, and this is precisely the aim of the 
present paper. This problem is also important in light 
of the recent progress in inducing superconductivity in 
graphene via the proximity effect, by building a Joseph- 
son junction,!^ as well as by growing a graphene sheet on 
rhenium .1^ 

The fact that superconductivity is not an intrinsic 
property of graphene has led to the search for the Dirac 
superconductor in materials with graphene-like proper¬ 
ties. A recently proposed type of engineered Dirac ma¬ 
terial consists of semiconducting nanocrystals with a 
truncated-cubic shape that self-assemble into a honey¬ 
comb superlattice.!^ The motivation to build these ma¬ 
terials was to study the electronic band structures that 
emerge when gapped, semiconducting systems are com¬ 
bined with features similar to those of graphene, by ar¬ 
ranging the nanocrystals in a honeycomb lattice. These 
materials have been experimentally realized for rock- 
salt PbSe and zinc-blende CdSe nanocrystals, which 
form honeycomb superlattices with a lattice parameter 
a 6nm via the attachment of the {100} facet of the 
nanocrystal.^ The electronic band structure has been 


theoretically described for PbSe, CdSe, and HgTe super¬ 
lattices, and it has been predicted that they exhibit Dirac 
cones in the conduction band above a wide gap enlarged 
by the quantum confinement 

Here, we investigate whether phonon-driven supercon¬ 
ductivity would be possible in artificial graphene samples 
made of PbSe or CdSe nanocrystals. We consider semi¬ 
conducting sheets that are either free-standing or capped 
with LiF, a dielectric which has been r ecentl y used to iso¬ 
late and passivate nanocrystal layers.^^^^ All of these 
systems are treated through an effective model where 
each nanocrystal is modelled as a superatom with a sin¬ 
gle effective s orbital representing the lowest cond uction 
state characterized by an s-envelope wavefunction.^^E^Hl 
We consider only the electrons close to the Dirac point. 
Furthermore, we assume that the electrons couple to a 
single Einstein phonon per superatom site, which corre¬ 
sponds to the longitudinal optical (LO) phonon for CdSe, 
PbSe, or LiF. This particular model is considered be¬ 
cause it can account for the features observed in scanning 
tunneling spectroscopy experiments on CdS^^or PbS^^ 
nanocrystals. 

In the problem studied in this paper, the LO phonon 
couples to the effective 5-electrons on the same site 
as well as on nearest-neighbor sites. The electrons 
are described by a tight-binding Hamiltonian, where 
the electron-phonon coupling includes both on-site and 
nearest-neighbor terms. We then integrate out the 
phonons to derive the effective electron-electron interac¬ 
tion. An estimate for this effective interaction is obtained 
based on a numerical analysis, and we find that for both 
PbSe and CdSe the effective interaction is repulsive, but 
can become attractive when the superlattice is capped by 
LiF. The renormalized values for the attractive on-site 
and nearest-neighbor interactions that we find indicate 
that in these materials only on-site pairing should oc¬ 
cur. Nevertheless, we theoretically investigate the more 
intricate case when nearest-neighbor interactions dom¬ 
inate over the on-site one, with the aim of motivating 
further experimental search for graphene-like materials 
that could fulfill these conditions and exhibit the elusive 
Kekule or hidden superconducting order described here. 

We consider both on-site and nearest-neighbor pair¬ 
ings, and the electron-electron interaction is then decou¬ 
pled in these channels using mean-field approximation. 
In this setup, we consider the problem of the compe¬ 
tition of the Kekule and the hidden supeconductors at 
both finite temperature and finite chemical potential in 
the vicinity of the Dirac points, and derive the gap equa¬ 
tions for these order parameters. First, excluding the 
possibility of the coexistence of the two orders, based 
on these equations, we find the dependence of the criti¬ 
cal (more precisely, crossover) temperature and the crit¬ 
ical interaction on the chemical potential, and analyti¬ 
cal solutions are obtained in the quantum-critical (strong 
coupling) and the BCS (weak coupling) limits, for both 
the Kekule and the hidden order. According to our re¬ 
sults, the Kekule order is preferred over the hidden order 


at both finite temperature and finite chemical potential. 
Second, when we allow for the possibility of coexistence, 
based on self-consistent mean-field gap equations, we ob¬ 
tain that above a critical value of the attractive interac¬ 
tion a mixed phase sets in where the two superconducting 
orders in fact do coexist. We show that the critical value 
for its onset is greater than the critical coupling for the 
hidden superconducting state in absence of the Kekule 
order. Therefore, there is a region in the phase diagram 
where the Kekule order is favored when competing with 
the hidden superconducting state. However, the latter 
eventually sets in and coexists with the Kekule super¬ 
conductor. Finally, according to our mean-field calcula¬ 
tions, the transition from the Kekule to the mixed phase 
is of the second order, but this result may change when 
fluctuations are taken into account. 

The paper is organized as follows. In Sec. [Hj we intro¬ 
duce the model that describes electrons, phonons, and 
their interaction in the system. In Sec. we define s- 
wave, Kekule and hidden superconducting order param¬ 
eters, and obtain the corresponding mean-field Hamil¬ 
tonians. We derive and analyze the gap equations for 
Kekule and hidden orders in Sec. hyi Results concern¬ 
ing the effective attractive interactions in self-assembled 
nano-crystals and the conclusions are presented in Sec.[V| 
Calculational details are presented in the Appendices. 

II. MODEL 

In this Section, we derive the effective model to study 
the superconducting properties of the system. We intro¬ 
duce full tight-binding and electron-phonon Hamiltoni¬ 
ans, after which we integrate out the phonon modes to 
obtain the effective model. 


A. Tight-binding model 

We describe the system using a tight-binding Hamil¬ 
tonian iLfuii, which includes nearest-neighbor hopping, 
Hubbard terms for electron-pairing on-site and between 
nearest-neighbors, the chemical potential, and electron- 
phonon coupling 

H'full = i^Hub + + H'el-ph. (1) 

The Hubbard Hamiltonian reads 

-f^Hub — ^ ^ ^ T h.C. 

+U y] + a^b 

i 

where (i,jf) denotes nearest-neighbor sites i and j, a] ^ 
{ai^cj) are creation (annihilation) operators for an elec¬ 
tron on sublattice A at site i with spin a, t is the 








nearest-neighbor hopping parameter, and U (V) are the 
on-site (nearest-neighbor) Coulomb interactions, respec¬ 
tively. The hopping parameter t depends on the size and 
the shape of the nanocrystals, as well as on the number 
of atoms connecting neighboring sites. T ypically, this pa¬ 
rameter is of the order of 10 The Hamiltonian 

for the chemical potential reads 


where /i < 0.51 to ensure the validity of the Dirac de¬ 
scription of the electrons. Finally, the electron-phonon 
Hamiltonian has the form 

^el-ph = ^ 

i 

^ H, (4) 

where the phonon frequency for phonons on sublattices A 
and B is equal, c]^ ■ {cA,i) create (annihilate) a phonon 
on sublattice A at site i, Sa connects nearest-neighbor 
sites with a = 1, 2, 3, and Vq and Vq are the coupling 
constants. 


where = 2n7r/{hP) for n G Z is the Matsubara fre¬ 
quency for bosons, 7k = 2N is the number of 

atoms in the system, and pA,q,cr,n is the electron density 

Pa, Cl, a,n — E i’A ,k+q, 

k,m 


The electron-phonon action for a phonon on sublattice 
B is obtained by the substitution A ^ B in the above 
equation. Completing the square and integrating out the 
phonon fields leads to the Hubbard terms U and V in 
Eq. @ renormalized by the electron-phonon coupling 

~ ^ ^ P^^n^\r^nPA, — ci,'\,—n A >■ B 

q,n 

-EE ^ PA,ci,a,nPB , — ci,a',—n’) (b) 

q,n a,a' 


where 
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The details of the derivation are presented in App. 
The effective Hubbard U (q) and V (q) are defined with a 
minus sign in the prefactor implying that for an attractive 
interaction these terms are positive. 


B. Effective model 

We now integrate out the phonons in Eq. Q, which 
results in additional contributions to the Hubbard terms 
U and V in Eq. ([^. To this end, we use 

2 : = jv[A,i’] J 


III. MEAN-FIELD APPROXIMATION 

Motivated by the possibility that the effective Hub¬ 
bard interactions may turn out to be attractive, we con¬ 
sider the superconducting instabilities in the artificial 
graphene samples. We first define superconducting order 
parameters, and then use the mean-field approximation 
to decouple the electron-electron interaction in Eq. §• 


where 0 {'ip) is an phonon (electron) field, the inverse Order Parameters 

temperature P = (/c^T)“^, and the action is given by 


S';(/)■*■, 0] = / dr[ilj^{r)dr'ip{r) + (l)^{r)dr(l){r) 

Jo 

+H ( 7 , </>)]• 

The electron-phonon action for the phonon on sublat¬ 
tice A yields 

^el—ph,A 0"^: 0 ] ~ ^ ^ ^A,q,n T 4^A,ci,r\ 

q,n 

T—^=== ^ ^ {^0 pA,ci,(T,n d- W 7q PB,q,cr,n^ 

^ (^A,-q,-n + 5 


The Hubbard Hamiltonian reads 

i 

(^) 

ihj) 


where it is assumed from now on that both IJ and V are 
positive. The order parameters corresponding to the on¬ 
site and nearest-neighbor pairing, respectively, have the 
following form 


( 5 ) 


Acr'^cr (^i:^j) — {bj,a'^i,a) 1 


( 10 ) 

( 11 ) 













where Aq represents the standard 5-wave order param¬ 
eter, and a general form for the nearest-neighbor order 
parameter Acr',cr is assumedP 

The electron densities in Eq. can be decoupled via 
a mean-field approximation 


ai,4,ai,t = Ao + , 

= Act',<t (rj, Tj) + 6 , 

such that 

H,s,u,v = 2Uj2\^of + vYl (i-i’ rOl' 

i {i,3) 

-u ^ (ai4ai,t + + ^.c. 

i 

—V (r^, Vj) bj^(j'(^i,a + h.c.^ , ( 12 ) 

(id) 


where qi = {qy^—Qx) and vf = 3at/2 is the Fermi veloc¬ 
ity. 

The Hamiltonian with the chemical potential in Eq. § 
isF^ = (l/2) with 

= -/i T3 0 (To 01. (15) 


2. On-Site Pairing Hamiltonian 


Only considering the terms in Eq. (12) that include the 
on-site order parameter Aq leads to = 4A[/|Aop + 
(l/2)Eq«'^MAo4' with 

Mao = -U [Re (Aq) Ti - Im (Aq) T 2 ] 0 ctq 0 ^ 7073- (16) 


3. Nearest-Neighbor Pairing Hamiltonian 


where the term quadratic in fluctuations O ((5^) is ne¬ 
glected. 


For the nearest-neighbor coupling, we use the Kekule 
ansatJ^ 


B. Hamiltonian in Dirac-Nambu Representation 


Aa,<T (r,, Tj) = A<^ cos (K ■ (r^ + r^)), (17) 


We now transform the full Hamiltonian in Eq. 0 with 
the U and V terms replaced by Eq. (12) to reciprocal 


space, expand around the Dirac points ±K, and use 
Dirac-Nambu representation to write the total Hamill- 
tonain as 


H = Eo + 


(13) 


where Eq is the energy of the condensate and the 16- 
component Dirac-Nambu spinors T’*' = with 

and arP 

“ (®K+q,<T ^K+q,<T ®T+q,<T ^T+q.tr) ’ 

^/i,cr(^) ~ (^K —q,cr ^K —q,cr K —q,cr ^ —K —q,cr) • 

The matrix M is given in terms of the 16 x 16 matrices 


^ijk = n 0cFj (g)7/c, 

where and aj are Pauli matrices acting in the particle- 
hole and spin space, respectively, and jk are 4x4 matrices 
acting in the sublattice-valley space defined as 70 = ctq ^ 
O's, 7i = OS (72, 72 = oq ^ (Ji, 73 = (7i (g) (72 and 
= 02 0 (72. 


1. Dirac Hamiltonian 


I (A;.t (r,, rj) + At4 (r^, rj))=A cos (K • (r^ + r^)), 

(18) 

) (li. ij) - (rt- rj)) = 11'. (19) 

where Eqs. 0 and represent a spin-triplet Kekule 
order, while Eq. ([T^ represents a spin singlet, the so- 
called hidden orderP In its full generality, the Kekule 
ansatz contains a phase. However, in the Dirac approx¬ 
imation we use here, i.e. only including electrons close 
to the Dirac points, the results are independent of this 
phase, and we have set it to zero. This degeneracy is, 
however, (weakly) broken when the lattice is reintro¬ 
duced in the problem.l^ 

Inserting this ansatz into Eq. ( p!^ leads to the following 
mean-field Hamiltonian for the Kekule order 

Tfxekuie = QNVm^ + i ^ (20) 

q 


with 

Mm = -V [{Xti - Ft2) 0 (73 + (I-T 2 - R-TQ 0 (Ti 
+ (Z+Ti + R+T2) 0 ( 72 ] 0 70- (21) 

Here, A = X + lY, R± = \ [Re (At) ± Re (A;)], J± = 
I [Im (At) ± Im (At)] and 

=X2 + y2+R^+j2 +ij^+/2_ (22) 


The hopping term in Eq. 0. can be written as Hd = For the hidden order, we find 
(1/2) with . 

Fa'= 1277F|Ay +(23) 


Md = VfTo®(Jo® «7o7i9i 





with 


Ma' = —V [Re (A') T 2 + Im (A') ti] ® (Ti ® i'yolsMD- 

(24) 

The proportionality to the Dirac Hamiltonian indicates 
that instead of opening a superconducting gap, the hid¬ 
den order renormalizes the Fermi velocity vp- 


IV. COMPETITION BETWEEN KEKULE AND 
HIDDEN ORDERS 

From now on, we only consider the Kekule and hidden 
order by setting the on-site order parameter to zero, i.e. 
Ao = 0. This case is more interesting to study because 
when all three superconducting orders are included, the 
s-wave order parameter is preferred, as is discussed in 
App. 0 We start by first considering the Kekule and 
hidden order parameters separately (excluding the pos¬ 
sibility of their coexictence) in the gap equations at both 
finite and zero temperature. Finally, we solve the self- 
consistent mean-field gap equations analytically at zero 
temperature and half filling. 


A. Gap equations 

First, we derive the thermodynamical potential for our 
system, which is followed by the computation of the gap 
equations and critical couplings. 


1. Thermodynamical Potential 


The thermodynamical potential ^ is obtained using 
the partition function 


z = = Tr 


(25) 


where iJtot is the sum of the Hamiltonians in Eqs. ( |14| )- 
(16), (20) and (23). Performing the trace yields 


^ = E 1^(1 




'0 


(26) 


q;s,s'=± 


where 


£;o = 6w(m2 + 2|Ay), (27) 


and Cjs^s' = sCjs' is obtained by diagonalizing + 

Mm + Ma' , with 


cos' = \j{vF\q\ + s'lif + |Ay, (28) 

where s = ±1 and s' = ±1 correspond to the spin and 
particle-hole degree of freedom, respectively. From this 
dispersion, we can see that the Kekule order m acts as a 




FIG. 1: (Color online) The dispersion relation in Eq. (28) is 
plotted for /a/t = A'/t — V/t — vp/t — \ and m — \ 
(panel (a)), showing that the Kekule order m opens a gap. 
In panel (b), the renormalization of the Fermi velocity by the 
hidden order parameter, given by Eq. (29), is shown. 


mass term for the Dirac fermions, and opens a supercon¬ 
ducting gap as shown in Fig.[^(a). The hidden order, on 
the other hand, for /i = 0, renormalizes Fermi velocity 
according to 


vf = tyj 1 + 


(29) 


as is displayed in Fig. (b). In Fig. the evolution of 
the thermodynamical potential with respect to different 
couplings is shown for both the Kekule and hidden order. 
We clearly see a second-order phase transition, with the 
thermodynamic potential at the critical coupling for the 
quantum phase transition shown in black. 











2. Finite-Temperature Gap Equations 


3. Critieal Coupling 


Minimizing Eq. (26) with respect to the Kekule order 


parameter leads to the finite-temperature gap equation 




s / 

S=±-^ 


'o — tanh —— 
(27r)2 uJa V 2 


V , (30) 


where cOg is given by Eq. (28). The equation determining 
the hidden order parameter at finite-temperature, which 
we loosely also call gap equation, is analogously obtained 


£4 

6 



dq \q? 

(27r)2 Cjg 


tanh 



(31) 



FIG. 2: (Color online) Evolution of the thermodynamical po¬ 
tential in Eq. ( [2^ for the Kekule (a) and hidden order (b) 
with = A = A/" = 1, ksTlt = 0.1, pi/t — 0. (a) Kekule 
order: V/t — 27 t (solid), V/t = Stt (black), V/t = Stt (large 
dashed), Vjt— IOtt (small dashed), and Vjt — 257r (dotted), 
(b) Hidden order: Vjt — IStt (solid), Vjt — IStt ( black), 
Vjt — 2l7r (large dashed), Vjt — SOtt (small dashed), and 
Vjt — SOtt (dotted). 


At half-filling (/r = 0) the Eermi energy is pinned at 
the Dirac points, and the density of states vanishes im¬ 
plying that there is a critical interaction at which the 
superconducting order sets in. Eor the Kekule order, we 
obtain the critical coupling 

-1 

, (32) 


y{T) = i-KVF 


A- 


21n(2) 

^Vf 


such that the zero-temperature critical interaction 
T4(0) = Vc = Sttvf/^- Eor the hidden order, we find 


VpT) = 


ISttF 

Vf 




3 ^<P) 


1 -1 


Fv 


(33) 


yielding for the zero-temperature critical interaction 
Vc{0) = = 187rt^/(i;i?A^), with the corresponding in¬ 

tegral for the hidden order solved in App. Here, A is 
the high-energy cutoff, up to which the continuum Dirac 
theory is valid, and which scales with the band-width of 
the order of the nearest-neighbor hopping. 

Taking VF/t = A = 1 and T = 0, we see that the criti¬ 
cal coupling for the Kekule order is smaller than the one 
for the hidden order. This implies that the system first 
enters the Kekule superconducting order and this state, 
therefore, dominates over the hidden order. In Eig.[^ we 
show how the critical interaction behaves with increasing 
temperature. We see that the critical interaction for the 
Kekule order remains smaller than that for the hidden 
order for any temperature. Eurt her more, the function 
Vc{T) {V^{T)) determines finite-temperature crossover 
from the quantum-critical semimetal to the Kekule (the 
hidden) superconducting state.l^Einally, we observe that 



FIG. 3: (Golor online.) Behavior of the critical interaction 
for the Kekule (blue) and hidden order (red) with increasing 
temperature with = A = 1, /a/t = 0, m/t = A' /t = 0. 





















the critical interaction for both orders increases with tern- At zero chemical potential, we find for the Kekule gap 
perature, since thermal fluctuations are expected to be 
detrimental for an ordered phase, meaning that the crit¬ 
ical coupling should increase. 


B. Zero Temperature 


The zero-temperature gap equations for the Kekule 
and hidden order are obtained from Eqs. (30) and (31), 
respectively, by setting T = 0 such that 


1 

1 


V ^ f dq 1 

3 ^^7 

V Vp f dq |gp 

6 (27r)2 Wj , ■ 


(34) 

(35) 


m (0, /u = 0) 


37rt>F A _ y \ 
yy v)' 


(37) 


Solving these integrals for weak and strong couplings, we 
can derive the zero-temperature gaps in both these limits. 


1 . Kekule Order 


To find the zero-temperature gap for the Kekule order, 
we set A' = 0 in Eq. (34), such that it simplifies to 


1 = 


V 






-l-/iln I - / = 

)/ 


—fi + Y /i^ -h K^m^(O)^ 


(36) 


Next, we solve the zero-temperature gap equation at 
finite chemical potential. Analytical solutions can only 
be found in the strong- and weak-coupling limit, V > Vc 
with m(0)/ju ^ 1 and V < Vc with m(0)/ju <C 1, respec¬ 
tively. Note that in the strong-coupling limit ju 1, so 
that the Eermi level is in the vicinity of the Dirac points. 
Therefore, this limit is governed by the quantum-critical 
point and we expect the resulting zero-temperature gap 
to exhibit power-law behavior. In the weak-coupling 
limit, on the other hand, the zero-temperature gap is 
expected to have a BCS-like form, since the system is 
away from the quantum-critical regime with a finite den¬ 
sity of states at the Eermi level. Applying these limits to 
Eq. (36) yields 


m (0, ju) 


m(0,^t=0) 

2 


1 + 



4/2^ 

mi0,fi=0f \ ’ 


^ exp[^m{0,fi = 0) - 1], 


V > y, m(0)//u > 1, 

V < Vc, m{0)/ii < 1. 


As expected, the superconducting gap shows quantum-critical power-law behavior at the strong-coupling, and is 
BCS-like in the weak-coupling limit. 


2. Hidden Order 


We 


solve Eq. (35) 


for m = 0 and obtain 


V vp I 

6 2tt 



A^ 

3Vl + 


jj? f 3 — \ 

4 


4a2 - 11 ( vTTc? 1 

-o -I- In - 

9 {l + a^) \ a 


2 - 3a2 \ 



where a = K|A'(0)|/t. Details of the calculation are zero chemical potential, we find 
presented in App. Eor the hidden order parameter at 


t 


1 - 



I A' (0,M = 0)| 


2 


(39) 
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FIG. 4: (Color online) Solutions of the finite-temperature gap equation for the Kekule (a) and hidden order (b) with vpjt — 
A = 1. In panel (a) IS! jt — 0 and V jt — dvr, while in panel (b) mit = 0 and V jt — 207r. In the plots we use the following 
values of the chemical potential: /i/t = 0 (solid), /i/t = 0.1 (large dashed), /i/t = 0.25 (small dashed), and [iji — 0.5 (dotted). 


which scales with t suggesting that this order parameter 
in fact does not open a gap, but renormalizes the Fermi 
velocity. 

At finite chemical potential, we find a solution in the 
strong- and weak-coupling limit, where |A'(0)|/t 1 

and |A'(0)|/t ^ 1, respectively, such that 


withF(yK') = ^(i-¥)-fg- 

C. Critical Temperature 


increases with increasing chemical potential, as expected 
from the fact that the density of states scales linearly 
with the chemical potential. An explicit expression for 
the critical temperature can be derived in the strong- and 
weak-coupling limit 

V > K, < 1, 

y < K, Pcii » 1, 

which again shows power-law behavior in the strong¬ 
coupling limit and BCS-like behavior in the weak- 
coupling limit. 


|A'|(0,a.) 



c>y, |A'(o)iA»i,, 
c<y, |A'(o)iA«i, 




2”' _ y } 21n(2)/cB 


+ /i 


Lastly, we determine the crossover (loosely called “crit¬ 
ical” hereafter) temperatures in the strong- and weak- 
coupling limit from the corresponding finite-temperature 
gap equations ( [M| ) and respectively, by requiring 

that m{Tc) = 0 (A'(T^) = 0) at the transition into the 
Kekule (hidden) order. 


1. Kekule order 


The solution of Eq. Pt at the critical temperature for 
A' = 0 is shown in Fig. |4(a)~ The critical temperature 


2. Hidden Order 


We solve the finite-temperature gap equation for the 
hidden order in Eq. ( [M] ) in a similar fashion by setting 
m = 0 and requiring ^( T/) = 0 at the transition. The 
solution is shown in Eig. |4(b)| with features similar to 
Eig. 4(a) The integral is solved in a similar fashion 
(outlined in App. and the explicit expression for the 
critical temperature reads 
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V > VJ, < 1, 

V < V', PciJi > 1. 


Fig. shows that the critical temperature for both the Kekule and hidden order parameters increases with chemical 
potential. Moreover, in Fig. |^the behavior of the critical temperature for both the Kekule and hidden order as a 
function of the coupling for different values of the chemical potential is displayed. We observe that the critical coupling 
decreases with increasing chemical potential. These features are expected on physical grounds, since the density of 
states linearly increases with energy as one moves away fro m the Dirac points. _ 



FIG. 5: (Color online) Behavior of the critical temperature 
for the Kekule (blue) and hidden order (red) with increasing 
chemical potential with vp/t — A = 1. For the Kekule order 
(blue), we chose V/t = Stt, and for the hidden order (red) 
V jt — IStt. 

D. Self-consistent gap equations at zero 
temperature and at half filling 

We now consider the competition of the Kekule and the 
hidden superconducting orders within the framework of 
the self-consistent mean-field gap equations (30) and (31) 


at zero temperature and zero chemical potential, which, 
after integrating over the angle, are rewritten as 


1 = 


Y_ 

Gtt 


I 

I 


1=1 “ 

Stt ./n \/q^ ^ Tn‘^ + A'p 


dq 


+ m? + A'p 


(40) 

(41) 


Here, we have conveniently redefined V^m? m? and 
v\V ‘^By rescaling the momentum, q 
g(l + |A'p)^/^, and performing the integration, we obtain 


V 


(42) 


l = ^(l + |Ar)-^(l-m), 

1 = T(l + |Ay)-2(l-2m2 + 2m3), (43) 

^ n 


where we redefined m/A ^ m. By inserting Eq. (42) 
into Eq. (43), and solving for the hidden order parameter 
we obtain six solutions of which only one is physically 
relevant. This solution has the form 


A'(yy,y) = 


V 

to 

(gi/s - 2VV^^ 

1 - GVsy]?' + ^2^2 + + Q2/3 

6143^1/3 


- 1 , 


(44) 


























FIG. 6: (Color online) Solutions of the finite-temperature gap equation for the Kekule (a) and hidden order (b) with vrjt — 
A = 1 for different values of the chemical potential: [iji — ^ (solid), /x/t = 0.1 (large dashed), /x/t = 0.25 (small dashed) and 
/x/t = 0.5 (dotted). m{Tc) = 0 and A'/t = 0 for (a), and A'(Tc) = 0 and m/t = 0 for (b). 


and yields the following Kekule gap 


m{V,V,,V') 


2Vc^ (wvi + Gi/ 3 ) + G^/^vv^ - _ 4^/ _ (52/3 


61 / 2 ^ 1/3 


(45) 


with the function G = G{V^Vc^V^) defined as 


G(V, Vc, vj) = + 4014 ® + + 33r2f^2-{7/2 _ 3y3f^/3 + 57 f^ 6 . 

,_ (461 



FIG. 7: (Color online) Behavior of the Kekule (blue) and 
hidden (red) orde r par ame ters with interaction strength V, 
given by Eq. (4^ and ( [44| , respectively. We use the critical 
coupling for the Kekule and hidden order in the absence of 
the competition W = 1 and K' = 5, respectively. 

As a result, above a critical value of the nearest- 
neighbor attraction Vc , we obtain a phase in which the 
Kekule and the hidden superconducting orders coexist. In 
Fig.[7| we plot the Kekule gap and the hidden order pa¬ 
rameter as a function of the nearest-neighbor attraction 
for fixed values of the critical couplings for the ’’bare” 
Kekule and the hidden orders, i.e, the critical couplings 


obtained without taking into account their competition. 
We observe that the value of the critical interaction Vc is 
greater than the bare value for the hidden superconduct¬ 
ing state. Hence, in the region V^ <V <Vc sketched in 
Fig.[^ the Kekule superconductor is favored over the hid¬ 
den order. However, the latter eventually sets in, and the 
two orders coexist. The dependence of the critical con- 











pling for the mixed phase on the bare critical coupling 
for the hidden order is shown in Fig. [91 We observe that 
Vc > expected based on the fact that the Kekule 
superconductor is gapped, while the hidden is gapless, 
which makes the former favorable over the latter. The 
preference of the Kekule order over the hidden order is 
so strong that for larger couplings, the system favors a 
mixed phase over a phase with a purely hidden super¬ 
conducting order. Furthermore, the transition from the 
Kekule into the mixed phase is of the second order, which 
consists of the two second order transitions in the sep¬ 
arate Kekule and hidden order channels. This feature 
may be an artifact of the mean-field approximation and 
when the fluctuations are included, this transition may 
turn out to be of the first order. However, this problem 
is beyond the scope of this work, and will be addressed 
in the future. 


V. DISCUSSION AND CONCLUSIONS 

Here, we show explicit values for the effective Hub¬ 
bard terms in Eqs. Q and (j^ obtained by describing 
the coupling of the electrons with the LO lattice defor¬ 
mations using a continuum dielectric model. The semi¬ 
conductor sheet is defined by its static 5in(0) and optical 
(high frequency) 5in(oo) dielectric constants. This ap¬ 
proximation i s com monly used to describe polarons in 
ionic material j^^ l ^^ l and L0- phonon coupling in semicon¬ 
ductor nano cry stals.^^^^^^^ The capping dielectric layer, 
if present, is described by 5out(0) and 5out(c>o). Details 
of the calculations are given in App. [E| 

Table jl] summarizes the results of the numerical calcu¬ 
lations performed on the superlattices described in Ref.|9l 
We consider nanocrystals with a truncated cubic shape 
and a size of 4.3 nm assembled in a honeycomb lattice. 
The (111) axis of the atomic lattice is oriented perpen¬ 
dicularly to the plane of the honeycomb sheet. The 
Coulomb interactions are obtained using 5in(0) = 10 and 
Sin(oo) = 6 for CdSe, and 5in(0) = 280 and ein(oo) = 25.2 
for PbSe. When the superlattices are free-standing, the 
effective interactions U and V remain repulsive, even 
though their magnitude is strongly reduced due to the 


Semimetal Kekule SC Mixed SCs 

-1 - 1- . I - > 

Vc K K” ^ 

FIG. 8: (Color online.) Phase diagram of the system as a 
function of the nearest-neighbor attraction V. As this cou¬ 
pling increases, at a critical value Vc the system first enters 
the Kekule superconducting state (blue region). In the re¬ 
gion labeled by dashed blue lines, the Kekule order is favored 
over the hidden order. The latter would in the absence of 
the Kekule order set in for V > VJ, but eventually coexists 
with the Kekule superconductor above the critical interaction 
Vc^ > Vc (red solid line). 


TABLE I: Parameters (in meV) defining the effective inter¬ 
actions in the free-standing superlattices of CdSe or PbSe, 
and in the super lattice of CdSe capped with LiF. Huje is the 
energy of the LO phonon which gives the strongest coupling 
to the electrons. 


System 

hjOE 

U 

V 

Vo 

Vo 

hpNU 

h/3NV 

CdSe 

26 

496 

262 

36 

7 

-360 

-666 

PbSe 

17 

290 

208 

32 

12 

-6 

-343 

CdSe/LiF 

82 

148 

66 

53 

27 

78 

9 



FIG. 9: (Color online.) Dependence of the critical coupling 
for the onset of the mixed phase Vc on the critical coupling 
for the hidden order in absence of the competition, W. The 
critical coupling for the Kekule order U/A = 1. 


coupling to phonons. In the case of PbSe, the effect is 
stonger because it is characterized by higher dielectric 
constants than CdSe. In fact, the effective interactions 
always remain repulsive due to the leakage of a large part 
of the electric field into the vacuum surrounding the su¬ 
perlattice when an electron is placed in a nanocrystal, 
implying that the dielectric screening from the ionic po¬ 
larization is unable to overturn the initial repulsive in¬ 
teraction. The opposite situation occurs when the CdSe 
superlattice is placed at 0.5 nm from a semi-infinite LiF 
sample [5out(0) = 8.9 and 5out(oc) = 1.9]. In that case, 
the effective interaction becomes positive due to the elec¬ 
tric field that strongly penetrates the external dielectric, 
and the main coupling comes from its polarization. The 
contribution coming from the polarization of CdSe can 
be neglected in a first approximation and huoE can be 
identified with its value in LiF. However, even in this 
case the on-site attractive interaction remains dominant 
over the nearest-neighbor one (see Table I and notice that 
positive values of energy actually correspond to attrac¬ 
tive interactions because of an overall minus sign in the 















Hamiltonian). Despite that, we theoretically investigated 
the more exotic superconducting orders that may arise 
when the nearest-neighbor attractive interaction domi¬ 
nates over the on-site one, which may also be relevant in 
the context of the ultracold atom systems.!^ 

More specifically, we have investigated the problem of 
the competition between the Kekule and hidden super¬ 
conducting orders in self-assembled artificial nanocrys¬ 
tals of graphene, at both finite temperature and finite 
chemical potential, within the mean-field theory for Dirac 
electrons, first by excluding the possibility of their coex¬ 
istence. As a result of this analysis, we find that the 
Kekule order is preferred over the hidden order at both 
a finite temperature and a finite chemical potential. On 
the other hand, within the self-consistent mean-field ap¬ 
proximation, allowing the coexistence, we find that there 
is a region in the phase diagram where Kekule order is 
favored as a result of the competition with the hidden 
superconducting order, but the latter eventually sets in 
and coexists with the Kekule state. Fluctuations may 
play an important role here, and addressing this problem 
requires the use of sophistic ated field-theoretical renor¬ 
malization group techniquesEven though our cal¬ 


culations suggest that if attractive interaction dominates 
in a self-assembled nanocrystal, it will be of the on-site 
type, they also indicate that there may be circumstances, 
as for instance even stronger screening, in which this re¬ 
sult could be overturned, so that the nearest-neighbor 
attraction could take over. This will hopefully motivate 
further search for materials where this will be the case, 
and would therefore open up a possibility for the realiza¬ 
tion of the exotic superconducting states in Dirac mate¬ 
rials. 
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Appendix A: Effective Hubbard terms 


The electron-phonon action for a phonon on sublattice A is given in Eq. (§• To integrate out the phonons, the 
square needs to be completed, which yields 


q,n a,a' ^ ^ ^ 


{VjQ Pa,C l,a, n T ^(q) PB,ci,(T,n) 




—ioJn + 


(u-O Pa, —q,cr', —n T u( Pb, — ci,cj', — n) 


j- n ^ ^ ^ ^ ^ (^0 PA,ci,a,n T P5,q,cr,n) ("^0 PA,—ci,a', — n T Q) P^, — q,cr' , — n) T A iA B . 

np ^^ T iUjOe 


q,n a,a' 


Here, uq = Vo/y/N and u(q) = Ko7q/v^ with 7k = Plugging this expression into the partition function 

leads to 


where 




= yp [A A] j 'D = JTi [A A] 

tkOE 

t cr,cr' ^ 

hijjE 


z = 


+ (^e) 


[('^0PA,q,cr,n T iSl)PB, ci,(t, n) (^OPA, —q,cr', —n T ^(^)P5, —q,cr'n)] 


\{uqPB, ci,<j,n T '^(^)pA,q,cr,n) ('^OPB, —q,cr'n T '^(^)PA,—q,cr',—n)] • 


q,n a,a' + (^eY 

Eqs. and (§ yield effective electron-electron on-site and on nearest-neighbor interactions in the form 

U _ huoE 


(Al) 


pq) = 


h/3 \N + {hjjEY 

hwE 




[p + t>(q)?;(-q)] | , 
[itov(-q) + Uov{q)] 


}■ 


h/3 \N’^ Cul^ihujEY 

Using that at finite temperature the zero Matsubara mode is dominant, and that |q|a <C 1 yielding exp (iq • S a) ~ 1 
such that 7 q 3, we then obtain the results in Eqs. Q and (H). 














Appendix B: Critical Interaction for Hidden Order 


To obtain the critical interaction for hidden order parameter in Eq. (33), we use 

3 1 

2 


(-e-2«)-±Li3 (-e-2“)+C, 


f 

/ duu^t8inh.{u) — d- i^^ln |l + e 
with C as a constant and Lin{x) is the polylogarithm function of the order n 


(Bl) 


Appendix C: Solving the Hidden Order Gap Equation 

To solve the hidden order gap equation, we use the following 


1 = 


s/ 

s=± 


dq 
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This integral can now be solved using Eq. (Bl). 


Appendix D: Critical coupling for s-Wave Superconductor 


Eor the s-wave order parameter, we find the following thermodynamical potentiaP 
=4iVf/|Ao|'-l ^ In \l+e^p[-ps^J{vF\q\+ s'^f 


^0| 

q;s,s'=± *- ^ 

Minimizing with respect to the s-wave gap Aq leads to the following finite-temperature gap equation 




1 


tanh 


^Ao;s 


2 


(Dl) 


where 


WAo;« = ^j{VF\q\+S|lf + (D2) 

The finite-temperature gap equation corresponds to the one for the Kekule order up to a prefactor, which is due to the 
gapped nature of both order parameters. Therefore, by minimizing the above thermodynamic potential and setting 
T = 0, we obtain the critical interaction for the 5—wave Dirac superconductor 


Uc = 271 Vp 



21n(2) 

I3vf 


(D3) 


We see that this critical interaction is smaller than those for the Kekule and hidden order, showing that the s-wave 
superconducting order is preferred. Since both the 5—wave and Kekule order parameter open up a gap at the Dirac 
points, the zero-temperature gaps and critical temperature are of the same form, and only differ in the prefactors. 


Appendix E: Effective couplings 

The parameters that define the effective interactions in Eqs. ([^ and ^ are obtained numerically by calculating 
the electrostatic interactions between electrons placed on the superlattices. The bare on-site Coulomb interaction of 
the Hubbard Hamiltonian is given by 


f/ = f/(oo) = I dV<^^(rW(r), 


(El) 






















where Pa{^) is the charge density corresponding to one electron in the s state of a nanocrystal A, and 0^(r) is 
the induced potential calculated by solving the Poisson equation using the high-frequency values for the dielectric 
constants of the inner and outer materials. Similarly, the bare nearest-neighbor Coulomb interaction is given by 

V = V(oo) = J (E2) 

where Pb(i*) is the charge density of an electron placed on a nanocrystal neighbor of A. For reasons that we 
clarify below, we also calculate U (0) and V (0) using the static dielectric constants instead of the high-frequency ones, 
thereby including the polarization coming from the LO phonons. 

The electron-phonon coupling terms Vb and Vo can be derived by writing the energy of the system in different 
electrostatic configurations. If we put one electron on a site on sublattice A, the classical energy derived from our 
model Hamiltonian is 

B(Q) = (E3) 

where Q = {2uje) is the operator corresponding to the displacement of the ions in response to the 

presence of the electron on sublattice A. The minimum of E{Q) at Qo = —V^Vq! gives the relaxation energy 

of the system after injection of the electron in nanocrystal A, the so-called Franck-Condon energy dpc = Vq /{Jiwe)- 
This energy is also given by the difference [t/(oc) — 1/(0)] /2 in the self-energy of the electron on a site on sublattice 
A in absence and in presence of the ionic response. The factor 1/2 comes from the adiabatic build-up of the charge 
in nanocrystal A. We then find 


Eo = y^[C/M-C/(0)]. (E4) 

In order to calculate Vq, we consider the quantity 1/(0) — H(0) which can be seen as the energy required to move 
a test charge (electron) from site H to A, when there is already an electron in A which induces the response of the 
ions. The analogue of this quantity derived from our model Hamiltonian is 


U{oc) + 2Ho 



H(oo) + (Ho + ^o) 



(E5) 


from which we find after some algebra, 


Vo = Vo 


V{oo) - E(0) 
U{oo)-U{0)' 


(E6) 


The quantities //(O), //(oo), H(0), and H(oo) are calculated numerically using the electron wave-functions directly 
derived from the atomistic tight-binding of Refs. mm- The charges on each atom (Cd, Pb, Se) are approximated 
by point charges, from which we find the potentials by integ rating the Poisson equation. _ 
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